In The Claims 



What: is clairaed is: 

Claim 1 (canceled) 

Claim 2 (currently amended) : The system estimation method 
according to claim ^1 or 9 or 11 , wherein the processing 
section calculates the existence condition in accordance with 
a following expression: 

= + ^fffH. > «. • = « * ('^> 

Claim 3 (currently amended) : The system estimation method 
according to claim ^1 , wherein the processing section 
calculates the existence condition in accordance with a 
following expression : 

-gEi + fn}>0, i = 0,---,A: (18) 

here, 

Q = 1-7), ^i = Y^^K7c ^ = ^-^('^^^ 

where the forgetting factor p and the upper limit value Yf have 
a following relation: 

Q < p = 1 - x.(y f ) ^ 1/ where jAl t ) denotes a monotonically 
damping function of Yf to satisfy = 1 and x(oo) =0- 

Claims 4-6 (canceled) 
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Claim 7 (currently amended) : The oyotcm cQtimation method 
according to claim 1/ — wherein at the otcp of executing the 
hyper 11^ filter, 

the proGcooincf acction calculatco thc A system estimation 

method for making state estimation robust and optimizing a 
fo rgetting factor p simultaneously in an estimation algorithm, 
in which 

for a state space model expressed by following 

expressions : 

Xk+I = F kX K + Gic Wk 

y^ = H kXk + Vk 

Zk = HkXic 

here , 

Xk: a state vector or simply a state, 

Wjc : a system noise, 

V)^ : an observation noise, 

yj^: an observation signal, 

Z}c: an output signal, 

F)^: dynamics of a system, and 

Gk: a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 
disturbance including the system noise w ^ and the observation 
noise Vk and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value Yf/ and 

the system estimation method comprises: 

a step at which a processing section inputs the upper 

limit value Yf/ the observation signal yk as an input of a 
filter and a value including an observation matrix Hk from a 
storage section or an input section; 



37 



a step at which the processing section determines the 

forgetting factor p relevant to the state space model in 
accordance with the upper limit value yt) 

a step of executing a hyper Hg^ filter at which the 

processing section reads out an initial value or a value 
including the observation matrix at a time from the storage 
section and obtains a filter gain Ks,k by using the forgetting 
factor p and a gain matrix Kk and by following expressions (20) 
to (22), or, expression (20) and expressions which are deleted 
Ji"-^ and Jl in the expressions (21) and (22), : 



x,^ = ji:,(:,i)/fle^(i,i) , k,=^pHp'''K,r:p:')JirI 





^ 1 

^k^k\k'i 


0 





Q{k) = 





0 







(20) 
(21) 

(22) 



Where, 

Rt = RijiRl, M = 



i 

p2 



0 



0 p^-yf 



1 0 
0 -1 



Sk\k-1 = -^^fclfc-l-Sfclfc-l 



Rejc = Rk + CtSkik-iCl, Ck = 



Hk 



. Rcji = RU'fiRljt, *o|o = «o 



(23) 



0(k) denotes a J-unitary matrix, that is, satisfies 
0(k)J0H(k)'^ = J, J = (Jl ® I), I denotes a unit matrix, 
Kk(:/1) denotes a coliomn vector of a first column of 

the matrix Kv . ^ 

v^hcrcin J ^' ^ and Jj, can be deleted in the oxprcooionct — (21) 
and (22) , 
here, 

x'^kik: the estimated value of the state Xk at the time k using 

the observation signals yo to yk/ 

yki the observation signal, 

Fk: the dynamics of the system, 

Ks,k: the filter gain, 

Hk: the observation matrix. 
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S'^kik: corresponding to a covariance matrix of an error of x'^kjk/ 
0(k) : the J-unitary matrix, and 
Re,k: an auxiliary variable. 

a step at which the processing section stores an 

estimated value of the state Xk by the hyper Ha> filter into the 
storage section; 

a step at which the processing section calculates an 

existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hj or 
the observation matrix Hj and the filter gain K^^j, and 

a step at which the processing section sets the upper 

limit value to be small within a range where the existence 
condition is satisfied at each time and stores the value into 
the storage section, by decreasing the upper limit value Yf 
and repeating the step of executing the hyper Hco filter. 

Claim 8 (currently amended: The system estimation method 
according to claim 7, wherein the step of executing the hyper 

Hoo filter includes: 

a step at which the processing section calculates and 

Z^k+i|k^^^ by using the expression (22) ; 

a step at which the processing section calculates the 
filter gain Kg^k based on feh^an initial condition and of S^ kjk:^ 
and an initial condition of Ck , and the matrix gain by using 
the expression (21) ; 

a step at which the processing section updates a filter 

equation of the Hoo filter of the expression (20) ; and 

a step at which the processing section repeatedly 
executes the roopcctivc otopo the step of calculating by using 
the expression (20) , the step of calculating by using the 
expression (21) and, the step of updating while advancing the 
t ime k . 
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Claim 9 (currently amended) : The ayatcm cotimation method 
according to claim 1, — wherein at the otep of executing the 
hyper H ^e filter, 

the proGcaaing gection A system estimation method for 

making state estimation robust and optimizing a forgetting 
factor p simultaneously in an estimation algorithm, in which 

for a state space model expressed by following 

expressions : 
Xk-fi = Fi^ X k + GkW)c 
y^ = H y.Xk + Vk 
Z}c = HkXk 
here , 

Xk: a state vector or simply a state, 

W k: a system noise, 

Vk: an observation noise, 

yk: an observation signal, 

Zk: an output signal , 

Fk: dynamics of a system, and 

Gk: a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 
disturbance including the system noise v/k ^nd the observation 
noise Vk and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value jf, and 

the system estimation method comprises: 

a step at which a processing section inputs the upper 

limit value Yf/ the observation signal yk as an input of a 
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filter and a value including an observation matrix Hy, from a 
storage section or an input section; 

a step at which the processing section determines the 

forgetting factor p relevant to the state space model in 
accordance with the upper limit value Yf/' 

a step of executing a hyper Ho^ filter at which the 

processing section reads out an initial value or a value 
including the observation matrix H } ^ at a time from the storage 
section and obtains the a filter gain Ks,k by using the 
forgetting factor p and a gain matrix Kk and by following 
expressions : 



^k\k = *Jb-l|fc-l + ^s^kiVk ~ -H'jfc^Jfe-llfc-l) 



0 



0 



(61) 
(62) 

e{k) (63) 



here, 0(k) denotes an arbitrary J-unitary matrix, and Ck 
= Ck+i^ is established, where 

Re^^Rk^ CkS^k\k^iCl, = [ ] » = Rlk^i^h* *o|o = *o (23) 
here, 

x'^kik: the estimated value of the state Xk at the time k using 

the observation signals yo to yk/ 

yk: the observation signal, 

Ks,k: the filter gain, 

Hk: the observation matrix, 

0(k): the J-unitary matrix, and 

Re,k: an auxiliary variable. 

a step at which the processing section stores an 

estimated value of the state X k by the hyper Hco filter into the 
storage section; 



41 



a step at which the processing section calculates an 

existence condition based on the upper limit value Yf 
forgetting factor p by the obtained observation matrix Hj or 
the observation matrix H j and the filter gain Kg^j, and 

a step at which the processing section sets the upper 

limit value to be small within a range where the existence 
condition is satisfied at each time and stores the value into 
the storage section, by decreasing the upper limit value Yf 
and repeating the step of executing the hyper Hop filter. 

Claim 10 (currently amended) : The system estimation 
method according to claim 9, wherein the step of executing the 
hyper Hoo filter includes: 

a step at which the processing section calculates K\ 
based on an initial condition of Re,k+i/ Rr,k+i and L'^k-hi b y using 
the expression (63) ; 

a step at which the processing section calculates the 
filter gain ICe^ — Kg^k based on the initial condition and by 
using the expression (62) ; 

a step at which the processing section updates a filter 
equation of the Hoo filter of the expression (61) ; and 

a step at which the processing section repeatedly 
executes the rcopoctivc otcpo the step of calculating by using 
the expression (63) , the step of calculating by using the 
expression (62) , and, the step of updating while advancing the 
t ime k . 

Claim 11 (currently amended) : The oyatcm estimation 
method according to claim 1, — wherein at tho step of executing 
the hyper H« filter. 
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the proccGoing PGCtion A system estimation method for 

making state estimation robust and optimizing a forgetting 
factor p simultaneously in an estimation algorithm, in which 

for a state space model expressed by following 

expressions : 
X)c4-i = F kXj^ + Gk Wk 
y i, = Hj cX k -t- Vk 
zj, = HkXk 
here, 

xj^ : a state vector or simply a state, 

Wk: a system noise, 

v^: an observation noise, 

y^: an observation signal, 

2k : an output signal, 

Fk: dynamics of a system, and 

Gk: a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 
disturbance including the system noise w^ and the observation 
noise Vk and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value Yf, and 

the system estimation method comprises: 

a step at which a processing section inputs the upper 

limit value Yf/ the observation signal yk as an input of a 
filter and a value including an observation matrix Hk from a 
storage section or an input section; 

a step at which the processing section determines the 

forgetting factor p relevant to the state space model in 
accordance with the upper limit value Yf/' 

a step of executing a hyper Ho : , filter at which the 

processing section reads out an initial value or a value 
including the observation matrix Hk at a time from the storage 
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section and obtains the a filter gain Ks,k by using the 
forgetting factor p and a gain matrix K"k and by following 
expressions : 



Xk\k = *ik-i|Jk-i + Ks^kiVk - HkXk^i\k-i) 
K,^, = p^Kk{:A)/Re,k{hl) 







0 


0 







0 

Kk 



ReMl - ^e,k ~ Ck-\-lLkRr]i^k^k-^l 
RrMl = - I^k^k-^l^e^k^M^k 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



Where, 



, = (w*+, u(fe + l-yv)] = [u(fc+l)«fc], H, =(u(l),0,...,0] 

ile., =il, + C'.i7i|„Cr, fl, = [j _^], i;,K, = diag{/,^p^...,p''+2J^ P=l-X(7/) 



1 0 
0 0 
0 1 



, 1^0 = 0, Xojo = sETo . Kk-p ^Kf^ (31) 



vvhorcin the abovo oxprcoDiono can be arranged aloo with 
rcopoGt to PCfe inotcad of K '^.^tT- 
here, 

yic: the observation signal, 

Fk: the dynamics of the system, 

Hk: the observation matrix, 

x'^kik: the estimated value of the state Xk at the time k using 
the observation signals yo to yk, 

Ks,k: the filter gain, obtained from the gain matrix K'k/ and 
Re,kf li'k^ an auxiliary variable. 

a step at which the processing section stores an 

estimated value of the state X k by the hyper Hop filter into the 
storage section; 
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a step at which the processing section calculates an 

existence condition based on the upper limit value Yf and the 
forgetting factor p by the obtained observation matrix Hj or 
the observation matrix H j and the filter gain Ks,i, and 

a step at which the processing section sets the upper 

limit value to be small within a range where the existence 
condition is satisfied at each time and stores the value into 
the storage section, by decreasing the upper limit value Yf 
and repeating the step of executing the hyper Hop filter. 

Claim 12 (canceled) 

Claim 13 (currently amended) : The system estimation 
method according to claim wherein an estimated value z\|k 

of the output signal is obtained from the state estimated 
value x^k|k at the time k by a following expression: 

z\|k - HkX'*'k|k* 

Claim 14 (currently amended) : The system estimation 
method according to claim wherein the Hoo filter equation 

is applied to obtain the state estimated value x^k| k==|'h^[k1 , Ll □ 

a pseudo-echo is estimated by a following expression: 

dk = Y.hi[k]uk~i. A: = 0,1,2,.-. (34) 

t=0 

and 

an echo canceller is realized by canceling an actual echo 
by the obtained pseudo-echo - 

Claim 15 (currently amended) : A system estimation program 
for causing a computer to make state estimation robust and to 
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optimize a forgetting factor p simultaneously in an estimation 
algorithm, in which 

for a state space model expressed by following 

expressions : 
Xk^i = F yXk + GjcW}, 
yx = HicXk + Vk 
Zk = Hk.Xk 
here, 

Xk: a state vector or simply a state, 

Wk: a system noise, 

Vki an observation noise, 

yk: an observation signal, 

Zk: an output signal, 

Fk : dynamics of a system, and 

Gk: a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 
disturbance including the system noise and the observation 
noise v > , and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value Yf/ ^J^d 

for a otatc apace model cxprcaacd by following 

Gxprcaoiono ; 

Y h I V fe 

^ = Hk -Xfe 
here, 

3€fe-i — a otatc vector or simply a otatc, 

w^-: — a oyotcm noioc, 

v^t-: — an observation noioc, 

^1 — an oboorvation oignal, 

— an output oignal, 
Ffe-: — dynamics of a oyotcm, — et¥tdt 
Gfe-i — a drive matrix. 
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a maximum energy gain to a filtor error from a 

diaturbancc weighted with the forgetting coefficient p ao an 
evaluation criterion io ouppreaaed to be omaller than a term 
corrooponding to a proviouoly given upper limit value y #-7 — and 

the system estimation program causes the computer to 
execute : 

a step at which a processing section inputs the upper 

limit value Vf , the observation signal yj^ as an input of a 
filter and a value including an observation matrix Hk from a 
storage section or an input section; 

a step at which the processing section determines the 

f orgettin g fac tor p relevant to the state space model in 
accordance with the upper limit value Yf; 

a step of executing a hyper Hop filter at which the 

processing section reads out an initial value or a value 
including the observation matrix at a time from the storage 
section and obtains a filter gain Ks.k by using the forgetting 
factor p and a gain matrix K"k and by following expressions: 
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K,,k = P^Kk{:.A)IR.*{iA) 







0 


0 







p-^LkR;}L,C,^i 

Re,k^k+lLk 



0 

Kk 



L J 
Where, 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



J 



A *+i = K+i u(A + 1 - iV)] = [u{k + 1) ujfc], ^1 = 0, . . . , 0] 



0 

LO -P7/J 



1 0 
0 0 
0 1 



-1 0 
0 /r^ 



, ^^0 = 0, Xo|o = xo, Kk=p~'K^ (31) 



here, 

Yjc: the observation signal, 

F)^: the dynamics of the system, 

H}c: the observation matrix, 

x^kik: the estimated value of the state xj, at the time k using 
the observation signals yp to y^, 

Ks.k: the filter gain, obtained from the gain matrix K'y^, and 
Re,k/ L^k' an auxiliary variable. 

a step at which the processing section stores an 

estimated value of the state Xk by the hyper Hoc, filter into the 

storage section; 

a step at which the processing section calculates an 

existence condition based on the upper limit value V f and the 
forgetting factor p by the obtained observation matrix Hj or 
the observation matrix H j and the filter gain Ks^j, and 

a step at which the processing section sets the upper 

limit value to be small within a range where the existence 
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condition is satisfied at each time and stores the value into 
the storage section, by decreasing the upper limit value Yf 
and repeating the step of executing the hyper Hoc, filter. 

Q 3tcp at which a proccooing ocction inputo the upper 

limit value Y #-7 — the obocrvation aignal y^ ^ aa an input of a 
filter and a value including an obaorvation matrix 1-1^ from a 
storage ocction or an input — ocction; 

a otcp at which the proccooing ocction dcterminco the 

forgetting coefficient p relevant to the otatc opacc model in 
accordance with the upper limit value js r 

a otcp at which the proccooing ocction reado out an 

initial value or a value including the obocrvation matrix IIfe --afe' 
a time from the otoragc ocction and uoco the forgetting 
coefficient p to execute a hyper filter cxprcooed by a 
following exprcooion ; 

= Fj e-^^e ^k i | k 1 I K e^4yfe Hfe-£V^3e^k i |k in - 
here , 

^n eotimated value of a otatc X k at a time k uoing 

obocrvation oignalo y^ to y ^— 
F'i^'-t — dynamicQ of — the oyotem, — and 
4 ^3 j e-? — ^ — filter gain, 

a otcp at which the proccooing ocction otorco an obtained 

value relating to the hyper filter into the otoragc 
ocction ; 

a otcp at which the proccooing ocction calculateo an 

cxiotence condition baoed on the upper limit value and the 
forgetting coefficient p by the obtained obocrvation matrix H^ , 
or the obocrvation matrix Hj, and the filter gain K ^-f — asad 

a otcp at which the proccooing ocction octo the upper 

limit value to be omall within a range where the cxiotence 
condition io oatioficd at each time and otorco the value into 



49 



the storage ocction by docroaoing the upper limit value and 
repeating the otop of executing the hyper li^ filter. 



Claim 16 (currently amended) : A computer readable 
recording medium recording a system estimation program for 
causing a computer to make state estimation robust and to 
optimize a forgetting factor p simultaneously in an estimation 
algorithm, in which 

for a state space model expressed by following 

expressions : 

Xk+i = F kX k + GicW ^c 

y>c = Hk XK + vj^ 

Zk = HkXK 

here , 

X)c: a state vector or simply a state, 

vjy^: a system noise, 

Vjg: an observation noise, 

yy^: an observation signal, 

Zjc: an output signal, 

F)^: dynamics of a system, and 

Gk: a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 
disturbance including the system noise W]c and the observation 
noise v^ and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value Yf, and 

for a otatc opacc model C3qprcoaod by following 

cxpreoaiona ; 
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here, 

— a otatc vector or oimply a atatc, 
w^-: — a oyatcm noioc, 
Vk--: — an obacrvation noiac, 
y^^-: — an obocrvation signal, 

— an output aignal, 

— dynamioo of a oyotcm, — asid 
Gk-: — a drive matrix, 

a maximum energy gain to a filter error from a 

dioturbance weighted with the forgetting coefficient p aa an 
evaluation criterion io ouppreaocd to be amaller than a term 
corroaponding to a prcviouoly given upper limit value y .in — arHdr 

the computer readable recording medium recording the 
system estimation program causes the computer to execute: 

a step at which a processing section inputs the upper 

limit value j f , the observation signal y : ^ as an input of a 
filter and a value including an observation matrix from a 
storage section or an input section; 

a step at which the processing section determines the 

forgetting factor p relevant to the state space model in 
accordance with the upper limit value Yf *' 

a step of executing a hyper Hop filter at which the 

processing section reads out an initial value or a value 
including the observation matrix at a time from the storage 
section and obtains a filter gain Kg^^ by using the forgetting 
factor p and a gain matrix K\ and by following expressions; 
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Xk\k = Xk-i\k-i + KsjiiVk - HkXk-i\k-i) 
K,.fc = /t»JKfc(:,l)/iie^(l,l) 



Kk+i 




0 


0 







T 



0 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



Where , 



. ^r*+i = K+itx(A + l-iV)] = (ii(/c + l)tijt], Hi = [w(l),0,...,0] 

P 0 



1 0 
0 0 
0 1 



0 -n} 



-1 0 

0 p-'^ 



, i7l|o = ciiag{p^p^....p^^2}, />=l-x(7/) 
, /fo = 0, Xo|o = aco » Kk - p ^K^ (31) 



here/ 



y^: the observation signal, 

Fk: the dynamics of the system, 

Hk: the observation matrix, 

x^kik: the estimated value of the state xy, at the time k using 
the observation signals yp to yk/ 

Ks,k: the filter gain, obtained from the gain matrix K'^f and 
Re.k/ l^k^ an auxiliary variable, 

a step at which the processing section stores an 

estimated value of the state Xk by the hyper Hcx> filter into the 

storage section; 

a step at which the processing section calculates an 

existence condition based on the upper limit value Y f and the 

forgetting factor p by the obtained observation matrix Hj or 
the observation matrix Hj and the filter gain Ks^i, and 

a step at which the processing section sets the upper 

limit value to be small within a range where the existence 
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condition is satisfied at each time and stores the value into 
the storage section, by decreasing the upper limit value Y£ 

and repeating the step of executing the hyper Hop filter. a 

otcp at which a proccDaing ocction inputo the upper limit 
value y ^-, — the obocrvation oignal input of a filter and 

a value including an obacrvation matrix Hj^ from a atoragc 
section or an input — ocction; 

a otcp at which the proccooing ocction dctcrminca the 

forgetting coefficient p relevant — to the otatc space model — ^ 
accordance with the upper limit value y^ -r 

a ptep at which the proccooing ocction reado out an 

initial value or a value including the obocrvation matrix H^ — 
a time from the otoragc ocction and uoeo the forgetting 
coefficient p to execute a hyper filter exprcooed by a 
following expreooion : 

here, 

x^fefi^ ; an eotimated value of a otatc x ^ at a time k uoing 
obocrvation oignalo y^ to y^ -r 
F-yr. — dynamico of the oyotem, — and 
¥z^^-t — a filter gain, 

a otcp at which the proccooing ocction otoreo an obtained 

value relating to the hyper filter into the otoragc 

ocction; 

a otcp at which the proccooing section calculateo an 

cxiotcnce condition baoed on the upper limit value y# 
forgetting coefficient p by the obtained obocrvation matrix H ^, 
or the obocrvation matrix Hj, and the filter gain K^ -^ — a«:d 

a otcp at which the proccooing ocction oeto the upper 

limit value to be omall within a range whore the cxiatcncc 
condition io oatioficd at each time and otoreo the value into 
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tho Dtoragc acction by dccrcaoing the upper limit value y ^—et^S: 
repeating the otcp of executing the hypor H.^ filter. 

Claim 17 (currently amended) : A system estimation device 
for making state estimation robust and optimizing a forgetting 
factor p simultaneously in an estimation algorithm, in which 

for a state space model expressed by following 

expressions : 
Xk^i = FkXk + GkW k 
yk = HkXK 4- vk 
Zk = HkXk 
here , 

Xk: a state vector or simply a state, 

Wk: a system noise, 

Vk: an observation noise, 

yk: an observation signal, 

Z k: an output signal, 

Fk: dynamics of a system, and 

Gv: a drive matrix, 

as an evaluation criterion, a maximum value of an energy 

gain which indicates a ratio of a filter error to a 
disturbance including the system noise Wk and the observation 
noise Vk and is weighted with the forgetting factor p is 
suppressed to be smaller than a term corresponding to a 
previously given upper limit value Vf, and 

for a atato apace model cxprcoacd by following 

cxprcooiona : 

3^; ^ - F k3€k I Gfe Wfe 

= H^ 3 ^ \ Vfe 

-25:56 — ss= Hfe5^ 

here , 

3€k-^ — a otatc vector or oimply a atatc. 
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wj^-: — a oyatcm noioc, 

^^-t — an obocrvation noioc, 

Y^". — an obacrvation oignal, 

^-z — an output oignal, 

F?^-: — dynamicD of a oyotctn, — afid 

Gie-: — a drive matrix, 

a maximum energy gain to a filter error from a 

dioturbancc weighted with the forgetting coefficient p aa an 
evaluation criterion io oupprooood to be omallor than a term 
corrcoponding to a previouoly given upper limit value . y g-T- 
the system estimation device comprises: 

a processing section to execute the estimation algorithm; 

and 

a storage section to which reading and/or writing is 
performed by the processing section and which stores 
respective observed values, set values, and estimated values 
relevant to the state space model, 
further comprising : 

a means at which the processing section inputs the upper 

limit value Yf/ t:he observation signal y^ as an input of a 
filter and a value including an observation matrix from the 
storage section or an input section; 

a means at which the processing section determines the 

fo rge tting fac tor p relevant to the state space model in 
accordance with the upper limit value Ye ; 

a means of executing a hyper Hop filter at which the 

processing section reads out an initial value or a value 
including the observation matrix at a time from the storage 
section and obtains a filter gain Ks,k by using the forgetting 
factor p and a gain matrix K\ and by following expressions; 
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Xk\k = *fc-i|*-i + KsjciVk - HkXk-i\k-i) 
Ks,k = P^Kki:A)/Red^A) 







0 


0 







0 

Kk 



— 1 ^ r 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



where , 



Ci|!+J = 



■ = u(A + l-yV)] = [u(fc + l)«t), =(«(!), 0,...,0] 



1 0 
0 0 
0 1 



, Ko = 0, xo|o = » Kk = p 'K^ (31) 



here^ 

Yk: the observation signal, 

Fk: the dynamics of the system, 

Hk: the observation matrix, 

x^kik: the estimated value of the state Xk at the time k using 
the observation signals yp to yk/ 

Ks,k: the filter gain, obtained from the gain matrix K'k/ and 
Re,k/ ^^k* auxiliary variable, 

a means at which the processing section stores an 

estimated value of the state Xk by the hyper Hop filter into the 
storage section; 

a means at which the processing section calculates an 

existence condition based on the upper limit value Yf and the 

forgetting factor p by the obtained observation matrix Hi or 
the observation matrix Hj and the filter gain Ks^j, and 

a means at which the processing section sets the upper 

limit value to be small within a range where the existence 
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condition is satisfied at each time and stores the value into 
the storage section, by decreasing the upper limit value Yf 
and repeating the means of executing the hyper H q q filter. 

the procc3aing ocction — input a the upper limit value y #t- 

thc obocrvation oignal y. ^ ao an input of a filter and a value 
including an oboorvation matrix Ilj^ from a atoragc section or an 
input — acct ion, 

the proGGoaing ocction dctcrminco the forgetting 

coefficient p relevant to the state space model — in accordance 
with the upper limit value y^ rr 

the proGCOoing scGtion reads out an initial value or a 

value including tho observation matrix at a time from the 
storage section and uaco the forgetting coefficient p to 
cxGGutQ a hyper H^^ filter cxpreased by a following exproooion ! 

F |^-,^3e ^jc ipc i— H — K9^-4^« H^Ffe-^^e ^c lik I n- 
here , 

— an estimated value of a state at a time 3c using 
observation signals y o to y ^— 

— dy nam ics of the system, — ajid 
Kp^k-; — a — filter gain, 

the processing section stores an obtained value relating 

to the hyper H, x > filter into the storage sec tio n, 

the processing section calculates an exi s tence condition 

based on the upper limit value y# ^^<^ t:hG forgetting 
coefficient p by the obtained observation matrix Hj, or the 
observation matrix Hj, and the filter gain Ke r^- — asd: 

the processing part sets the upper limit value to be 

small within a range where the existence condition is 
satisfied at each time and stores the value into the storage 
section by decreasing the upper limit value y# repeating 
the step of executing the hyper filter. 
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Claim 18 (new) : The system estimation method according to 
claim 9, wherein the processing section calculates the 
existence condition in accordance with a following expression: 

-gEi + /97j>0, i = Or-,k (18) 

here, 

Q = 1-7?, ^i = Y^^K7c ^ = (19) 

where the forgetting factor p and the upper limit value yf have 
a following relation: 

0 < p = 1 - x(Yf) - 1' where x(yf) denotes a monotonically 
damping function of yf to satisfy x{l) = 1 and x{oo) = 0. 

Claim 19 (new) : The system estimation method according to 
claim 9, wherein an estimated value z\|k of the output signal 
is obtained from the state estimated value x'^kjk at the time k 
by a following expression: 

z\|k - HkX'^klk- 

Claim 20 (new) : The system estimation method according to 
claim 9, wherein the Hoo filter equation is applied to obtain 
the state estimated value x\ik=[h^i[k],nDn, h^[k]] 

a pseudo-echo is estimated by a following expression: 

dk = $3/ii[A;K_i, A: = 0,l,2,-- (34) 

i=0 

and 

an echo canceller is realized by canceling an actual echo 
by the obtained pseudo-echo. 

Claim 21 (new) : The system estimation method according to 
claim 11, wherein the processing section calculates the 
existence condition in accordance with a following expression: 
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-gEi + P7/>0, i = 0,--,A: (18) 

here, 

Q = 1-7?, ^ = Y^^k7^ P=1-x(7/) (19) 

where the forgetting factor p and the upper limit value yf have 
a following relation: 

0 < p = 1 - X{Yf) ^ 1^ where x(Yf) denotes a monotonically 
damping function of Yf to satisfy %(!) = 1 and x{co) = 0. 

Claim 22 (new) : The system estimation method according to 
claim 11, wherein an estimated value z\|k of the output signal 
is obtained from the state estimated value x'^kik at the time k 
by a following expression: 

Claim 23 (new) : The system estimation method according to 
claim 11, wherein the Hoo filter equation is applied to obtain 
the state estimated value x\|k=[h'^i[k],DDD, h^[k]] 

a pseudo-echo is estimated by a following expression: 

N-l 

dk = E^^il^K-i' A: = 0,1,2,.-. (34) 

i=0 

and 

an echo canceller is realized by canceling an actual echo 
by the obtained pseudo-echo. 
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